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® Reliable noninvasive measurement of blood gases. 

@ This Invention relates to methods and apparatus for, preferably, determining non-invasively and in vivo at 
least two of the five blood gas parameters (i.e., pH, [HCOa'], PCO2, PO2, and O2 sat.) in a human. 
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BACKGROUND OF INVENTION 

This invention relates to both a method and apparatus (as illustrated In Figure 1) for, particularly, the 
noninvasive determination of those parameters currently reported on a standard clinical arterial blood gas 
5 report, namely: hydrogen ion concentration (pH) partial pressure of carbon dioxide (PCOa). partial pressure 
of oxygen (PO2), bicarbonate concentration ([HCO3"]) and oxygen saturation (O2 sat). 

Arterial blood gas determination is the cornerstone of diagnosis and management of cardiopulmonary 
disease in the critically ill patient. As effective oxygenation and maintenance of acid-base balance in such a 
patient is necessary for survival, measurement of arterial blood gases is typically the most frequently 
70 ordered laboratory test in a hospital's Intensive care unit. 

The standard arterial blood gas report contains the following information: pH, PCO2, PO2. [HC03~], and 
O2 saturation. The pH and PCO2 provide valuable information regarding acid-base and ventilation status. 
The bicarbonate level provides additional infonmatlon on acid-base balance which allows the physician to 
determine whether an acid-base abnormality is respiratory or metabolic in origin. The two other indices. 
75 PO2 and O2 saturation, reflect the amount of oxygen present in the patient's blood. 

At the present time standard clinical practice requires arterial puncture for procurement of an arterial 
blood sample. The arterial puncture is painful to the patient and associated with a variety of complications. 
Additionally, the process for analyzing the arterial blood sample is lengthy, requires multiple personnel and 
is not continuous. 

20 Due to the important clinical role of the information obtained by arterial blood gas analysis and the 
drawbacks of conventional clinical analysis, a number of alternative technologies have been proposed to 
measure one or more blood gas parameters. The prior art can be divided into 7 major categories. 

1. Invasive colorlmetric technologies requiring direct contact between the colorimetric substances and 
arterial blood. These technologies measure pH, PO2, and PCO2. 
25 2. Separation of the blood gases by semipermeable membranes, with subsequent concentration 
determination by absorption spectroscopy. The technology is limited to the measurement of PO2, and 
PCO2. 

3. Electrode devices for measurement of PO2, and PCO2 in peripheral skin. This technology has been 

commercially developed by Radiometer, Denmark. 
30 4. Invasive transistor devices for measurement of pH, and PCO2 in blood. 

5. Semi-invasive techniques for measurement of pH, PCO2 and PO2 in peripheral skin. 

6. Optical or spectroscopic determination of PCO2 and PCO (partial pressure of carbon monoxide) in 
tissue; 

7. Pulse oximeters for noninvasive measurement of arterial blood oxygen saturation (O2 sat.). 

35 The terminology used in clinical medicine, technical publications and prior patents is not always 
consistent. In the clinical practice of medicine the term "blood gas" is specific for the determination of 
certain factors in arterial blood. The determination of blood gases in venous blood is of little or no clinical 
significance and is not standard practice. Technically the term "blood gas" is misleading as not ail of the 
components listed on a standard laboratory blood gas report are specific for gas measurements. PH and 

40 [HCOa"] are concentration measurements of non-gaseous substances. In this application "blood gas 
parameters" are pH. PCO2. PO2, or [HCOa'] and O2 saturation. 

The foregoing parameters are not always measured in arterial blood; they can be measured in tissue 
and skin. Thus, it is important to distinguish between measurement in arterial blood, tissue and peripheral 
skin, as their sensitivity and clinical utility differ. Blood gas measurements in arterial blood are currently the 

45 reference standard used in clinical medicine. 

Several prior art methodologies measure certain blood gas parameters in the peripheral skin or tissue. 
Peripheral skin and tissue are differentiated in that skin is the outermost layer of the body, while tissue 
refers to the sum total of body mass existing between two skin surfaces. Tissue measurements, if accurate, 
may be the measurements of highest clinical utility as such measurements define the acid-base and oxygen 

50 status of the peripheral tissue. The maintenance of normal physiology in these tissues is the goal when 
treating cardiopulmonary disease. 

Measurement of blood gas parameters in peripheral skin has limited clinical utility. When peripheral skin 
measurements are compared to arterial blood measurements, it has been found that the skin measurements 
of PO2 are not adequate. Peripheral skin measurements are also less sensitive to acute changes in blood 

55 gas concentrations. 

To summarize, there are five blood gas parameters which can be determined in three different media, 
(1) arterial blood, (2) tissue, and (3) peripheral skin. A ranking of the current clinical utility of these 
measurements would be arterial blood, tissue and peripheral skin a distant third. 
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Invasive Coiorimetric Technologies 

The application of coiorimetric measurement techniques to the measurement of pH. PCO2 and PO2 is 
well summarized by: 

5 U.S. patent No. 4,854.321 titled "Integrated Optical System for Monitoring Blood Gases" by A. A. 
Boiarski describes the use of coiorimetric substances in dye wells. A coiorimetric substance is a chemical 
which changes color in a quantifiable way when exposed to another given substance. The technology 
requires that the probe containing the dye wells be in direct contact with the blood. 

U.S. patent No. 4,682.985 titled "Fiber Optic Probe for Quantification of Coiorimetric Reactions" by D. 

10 Costello discloses a similar technology with the additional feature of a probe which has an overall diameter 
"sufficiently small to permit the probe to be inserted into living tissue directly" (column 1, lines 12-17). A 
coiorimetric substance contained in the sample chamber changes colors in response to chemical properties 
of the chemical to be colorimetrically measured, thereby changing the amount of light transmitted through 
the sample chamber by the optical fibers. 

75 A similar technology to that disclosed by Boiarski and Costello Is described for pH detenmlnation by S. 
R. Goldstein et al. in "A Miniature Fiber Optic pH Sensor For Physiological Use." Journal of Biomechanical 
Engineering. Vol. 102, 141-146. The article describes the concept of using two fiber optic strands to 
illuminate and remotely sense the color change of a dye indicator contained in a hollow fiber permeable to 
hydrogen ions. 

20 Another coiorimetric technology is described by Gehrich et al. in U.S. patent No. 4, 989,606. The patent 
describes an invasive system for sampling small amounts of blood from a patient in a continuous manner. 

Separation by semipermeable membranes. 

25 The direct measurement of the partial pressure of a PO2 and PCO2 gas by absorption spectroscopy 
following separation by a semipermeable membrane is described in U.S. patent No. 4,201 ,222 to T, Hasse. 
The technology is based on the concept that gas can flow through a barrier until an equilibrium state is 
reached. If the barrier Is made of a material permeable to selected gases then absorption spectroscopy can 
be used to quantify the amount of gas in the chamber. The specific method of measurement uses a single 

30 wavelength determination. The device described requires insertion into an arterial blood vessel and is not 
suitable for noninvasive measurements. 

Eiectrode Devices/Sensors 

35 A generalized description and overview of this technology can be found in P. Rolfe's article "Review of 
Chemical Sensors for Physiological Measurement." Biomedical Engineering Centre, University of Oxford, 
UK. Journal of Biomedical Engineering. 1988, April, 10(2), pages 138-45. The application of this technology 
to transcutaneous or noninvasive blood gas determination Is most completely described in a series of 
patents assigned to Radiometer A/S of Denmark. U.S. Patent Nos. 4.274,418; 4,324,256; and 4,685.465 

40 describe the use of potentiometric and polarographic electrodes for measurement of the blood gas partial 
pressures of oxygen and carbon dioxide. It is important to note that electrodes measure parameters of the 
medium In direct contact with the probe. Thus, the partial pressures of oxygen and carbon dioxide 
measured by the transcutaneous electrodes would be that of the peripheral skin, not arterial blood. 

A similar electrode based instrument is described by G. J. Ullrich et al. in U.S. patent No. 4,930,506. 

45 The technology varies from the Radiometer patents in that the Ullrich invention couples an electrode device 
with a standard pulse oximeter for continuous detemnlnation of peripheral skin PCO2 and arterial blood 
oxygen saturation. 

In-Vivo pH Measurement Transistors 

50 

M. E. Meyerhoff et al., U.S. patent No. 4.694,834, B. Oeseburg et al. "Intravascular pH-ISFET, a method 
of the future," Scand. J. Clin. Lab Invest 1987 (47. Suppl. 188; 31-35). and P. Bergueld "The Development 
and Application of FET-based Biosensors," Biosensors 2 (1986) 15-33, describe the use of Ion selective 
field effect transistors (ISFET), specifically, potentiometric ion selective electrodes based on semiconductor 
65 technology for the in-vivo continuous determination of pH. The technology is currently limited to pH 
determination and is invasive. 
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Semi-Invasive Techniques 

A seml-lnvasive technology for the measurement of pH, PCO2 and PO2 is described by M. A. Fostick in 
U.S. patent No. 4,041,932. The technique includes the formation of a skin 'window' in a small area of a 

5 patient by removal of substantially all of the stratum corneum (top dry layer of skin) in the area. An 
enclosed chamber is positioned tightly against the skin around the window so formed and sealed thereto. 
The chamber and adjacent skin is heated to above normal skin temperature. Gases or fluids are collected In 
the chamber until equilibrium is reached. PH. PCO2, and PO2 are then measured by absorption spec- 
troscopy via univariate methodology. This technology determines these blood gas parameters in peripheral 

10 skin, not arterial blood. 

Optical or Spectroscopic Determinations 

U.S. patent No 4.509.522 by T. J. Manuccia describes the detemlnation of PCO2 and PCO; invasively 
75 in arterial blood, or non-invasively in tissue. The gas concentration is determined by a univariate linear 
algorithm. Univariate algorithms cannot work in a situation where there are varying analytes with overlapping 
absorptions. Although this methodology may work well for CO2 and CO, it would not work well for 
determination of pH. PO2, [HC03~] and O2 saturation. The spectral regions used for determination of pH, 
PCO2, and [HC03~] overlap extensively. Thus, the univariate method described by Manuccia et al. is not 
20 applicable to the simultaneous noninvasive determination of all blood gas parameters. 

The clinical utility of Manuccia's technology will be limited due to the frequencies used for the PCO2 
and PCO measurement. The frequencies specified by Manuccia (I.e. between 3.0 and 14.0 microns) are in 
the mid-infrared, an area of the spectral domain in which the light propagation characteristics are not 
consistent with transmission measurements through the finger, and the penetration of light is not to the 
26 arterial vasculature. Thus, measurements made will be similar to a peripheral skin determination. To 
summarize, Manuccia's technology cannot measure all blood gas parameters, uses frequencies inconsistent 
with noninvasive tissue determination and employs an algorithm which will perform poorly in the near 
infrared frequency region. 

30 Oximetry 

Pulse oximetry is an accepted method of oxygen determination and has been utilized in clinical 
medicine for years. All pulse oximeters are based on several governing principles. First, the concentration 
of blood In a given location of the body varies with each pulse of the heart. With each heart beat a systolic 

35 pulse pressure is generated which leads to a maximal expansion (i.e. dilation) of the vascular system. 
During the resting period of the cardiac cycle (i.e., diastole) there is no pressure generated and the vascular 
system returns to a minimal size. The light transmitted or reflected during diastole interacts with the skin, 
fat. bone, muscle and blood. Light transmitted or reflected during systole interacts with the same skin, fat, 
bone, muscle, and blood, plus an additional amount of blood which is present due to the expansion of the 

40 arterial system. If the diastolic signal is subtracted from the systolic signal the result is a signal which 
represents the additional amount of blood. The subtraction process removes the interferences created by 
the Interaction of the light with the skin, fat, bone and muscle. The quality and clarity of the subtraction 
generated signal Is related to the amount of additional blood present, which Is proportional to the pulse 
pressure. See Figure 2 for a graphical representation of the above process. 

45 All present pulse instruments assess variations in red blood cell concentration by utilizing a light 
frequency near or at the isosbestic point, where optical measurement of pulsatile volume is made 
Independent of oxygen saturation. An isosbestic wavelength is one which changes only with blood 
concentration but does not change intensity with oxygen saturation. Consequently, such a wavelength 
(typically in the range of 800-850 nm) Intentionally eliminates Information about oxygen saturation and 

50 establishes a reference. A second wavelength in the red portion of the spectrum, which is sensitive to 
oxygen saturation, is detected by either a transmission or reflection sampling technique. By using the 
isosbestic wavelength as a reference and by comparing Its spectral Intensity to the intensity of the second 
wavelength in the red portion of the spectrum, it is possible to determine the oxygen saturation of the 
arterial blood non-lnvasively. 

55 
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Spectroscopic Determination of pH in Non-Biological Systems 



Optical determination of pH in non-biological systems has been demonstrated by several investigators. 
Two representative articles are "Near-infrared Spectrometric Determination of Hydrogen Ion. Glucose and 

5 Human Serum Albumin in a Simulated Biological Matrix," by Drennen et al. and "Measurement of caustic 
and caustic brine solutions by spectroscopic detection of the hydrogen Ion In the near-infrared region, 700- 
1150 nm". by M. K. Phelan. 

With specific reference to pH determination, Drennen et al. describe a process for quantitative 
spectroscopic measurement using smoothed second-derivative spectra with multivariate analysis by princi- 

10 pie component regression (PGR). Although the paper does describe the quantitative spectroscopic deter- 
mination of pH by multivariate analysis, the solutions were extremely simple with the pH being the sole 
source of variation. This is sharp contrast to determination of pH in blood or tissue which are extremely 
complex systems with multiple varying parameters. (i.e., the remaining four blood gas parameters). Also, 
the pH range measured by Drennen et al. is outside the range consistent with human survival and the 

75 frequency region used is not compatible with noninvasive transmission measurement. Specifically, fre- 
quencies greater than 1450 nm are strongly attenuated by water and transmission through a significant 
amount of tissue is not possible. 

Phelan et al. demonstrated measurement of hydrogen ion concentration in caustic (sodium hydroxide) 
and caustic brine (sodium hydroxide and sodium chloride) solutions. These investigators used the fre- 

20 quency region from 700-1150 nm with subsequent analysis of the resulting spectra by multiple linear 
regression and partial least squares. As above, it's important to emphasize that the system studied is not 
biological, the pH variation Is far outside that of normal physiology, and is not nearly as complex as found in 
biological systems. 

25 Bacl(ground on Spectral Analysis Algorithms 

Optical measurements in whole blood following reactions with enzymes and/or reagents are commonly 
employed In the clinical chemistry laboratory. The use of spectroscopic techniques with multiple 
wavelengths for such a determination is much less common. Further, it is essential to realize that all prior 

30 noninvasive monitors, specifically Manuccia, et al.. and the majority of prior art pulse oximeters used 3 or 
less measured intensities and one or two variables for analysis. Methods that simultaneously use two or 
more variables are known as multivariate methods. Thus, while Manuccia et al and the above described 
pulse oximeters do utilize optical measurement techniques, the algorithms used for analysis are not as 
powerful or sophisticated as those utilized herein. 

35 A simple illustration of the increased capability of multivariate methods in component concentration 
determination Is provided by Figures 3. In Figure 3A one can see that an Impurity component, whose 
spectrum overlaps that of the analyte, can affect the spectrum of the analytic band. Therefore, the accuracy 
of the analysis will suffer when the analysis is performed at a single wavelength v\ or when rationing v^ to a 
reference wavelength. The measured absorbance, A^, at the analysis wavelength, fi, for a sample 

40 containing the impurity is different than the true absorbance. At, of the analyte at that wavelength. If the 
calibration curve in Figure 3B is from spectra of samples containing no impurity, then the presence of the 
impurity in the sample will yield an apparent concentration that may be quite different from the true 
concentration. This error will remain undetected if the intensity was measured at only one wavelength. If the 
impurity is included in the samples, a calibration plot similar to that in Figure 38 will exhibit large scatter 

45 among the data, and the result will be both a poor calibration curve and concentration estimates that have 
poor precision for the unknown samples. However, with analysis at more than one wavelength, not only can 
the presence of the impurity be detected, Figure 30, but if Its presence is included in the calibration, 
quantitative analysis of the analyte Is possible with multivariate calibration methods, even if the Impurity and 
its concentration are unknown. 

50 An indication that the unknown Is different from the set of calibration samples not containing the 
impurity is obtained by plotting the absorbance of the calibration samples and the unknown sample spectra 
at two frequencies selected for analysis. As exhibited in Figure 30, the spectrum of the sample containing 
the impurity (Indicated by "x") is obviously different than that of the calibration spectra (i.e., it Is an outlier). 
Outliers are those samples or spectra among either the calibration or unknown data that do not exhibit the 

55 characteristic relationship between composition and spectra of the other calibration samples. The sensitivity 
in detecting outliers is Increased by increasing the number of frequencies Included in the analysis. The 
number of independently varying Impurities that can be simultaneously accounted for in the analysis Is also 
increased by increasing the number of frequencies utilized. 
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Accurate univariate methods are dependent upon the ability to identify a unique, isolated band for each 
analyte. Multivariate methods can be used even when there is overlap of spectral information from various 
components over all measured spectral regions. Unlike univariate methods, multivariate techniques can 
achieve increased precision from redundant information on the spectra, can account for base-line variations, 

5 can more fully model nonlinearities, and can provide outlier detection. 

The general approach that is used when statistical multivariate methods are applied to quantitative 
spectroscopy problems requires calibration in which a mathematical model is generated relating analyte 
concentrations to reference spectra. See Figure 4, This calibration model can then be used for prediction of 
concentrations in unknown samples. The spectra of a series of calibration standards are first obtained, such 

10 that the spectra span the range of variation of all factors which can influence the spectra of future unknown 
samples. Assuming that the calibration uses samples that contain all the components expected in the 
unknown samples and spans their expected range of variation, the calibration will be able to empirically 
account for (or at least approximate) non-ideal behavior in Beer's law, Independent of the source of the non- 
ideal behavior. As used in this application "nonlinear" refers to any deviation in Beer's law or the inverse 

75 Beer's law relationship (i.e.. which cannot be modeled with the standard linear expression y = mx b; 
where y represents the dependent variable, x Is the independent variable, and m and b are, respectively, 
the slope and intercept). 

Once the empirical calibration relating spectra and component concentrations has been performed, then 
the spectrum of the unknown sample can be analyzed by a multivariate prediction step to estimate the 

20 component concentrations or properties. In addition, spectral residuals (i.e., the difference between mea- 
sured and estimated spectra) can be used to determine if the unknown sample spectrum is contained within 
the range spanned by the calibration samples. If the unknown sample is not representative of the calibration 
samples (i.e., Is an outlier), spectroscopic interpretation of the residuals can often be made to determine the 
source of any differences between unknown and calibration samples. See Haaland, David M.: "Multivariate 

25 Calibration Methods Applied to Quantitative FT-IR Analysis" In Practical Fourier Transform Infrared Spec- 
troscopy, Industrial and Laboratory Chemical Analysis, Edited by J.R. Ferraro and K. Krishman, Academic 
Press, Inc. 1990. Not only do multivariate statistical methods provide enhanced analysis of component 
concentrations, but such multivariate methods have also recently made possible the estimation of physical 
and chemical properties of materials from their spectra. 

30 The use of multivariate analysis In noninvasive medical monitoring is best described in U.S. patent No. 
4,975,581 to M. R. Robinson, K. J. Ward, R. P. Eaton and D. M. Haaland and current pending European 
patent application, Publication No. 0 522 674 A2, titled "Oximeter for Reliable Clinical Determination of 
Blood Oxygen Saturation in a Fetus" by M. R. Robinson. D. M. Haaland and K. J. Ward. U.S. patent No. 
4,975,581 discloses a method and apparatus for, particularly, quantitatively determining the amount of 

35 glucose in a human. 

A number of multivariate calibration methods are available for quantitative spectral analyses. Many of 
these have been reviewed by D. M. Haaland, supra. These include PLS, PCR. CLS, MLR, Q-matrix, Kalman 
filtering, and cross-correlation. In addition, ridge regression, continuum regression, and neural networks are 
other possible multivariate methods that can be used in quantitative spectral analysis. 

40 Multivariate methods which are well suited for analysis of blood gas parameter spectroscopic data are 
those that model the spectra using an inverse Beer's law model, such as principal component regression 
(PCR) or partial least squares (PLS). An advantage of this multivariate approach is that the nonlinearities in 
the spectral response to changes in composition can be accommodated without the need for an explicit 
model. PLS and PCR methods are capable of achieving accurate and precise results in the presence of 

45 linear and nonlinear dependencies in the absorbance spectrum at various frequencies. Thus, an entire 
spectral region can sometimes be used in multivariate analysis without the requirement that the spec- 
troscopist choose an optimal set of wavelengths for the analysis. Similariy, these methods of computation 
are not sensitive to linear dependencies introduced by over sampling of information at many frequencies in 
the construction of the calibration samples. See Cahn. et al., "Multivariate Calibration of Infrared Spectra for 

50 Quantitative Analysis Using Designed Experiments". Applied Spectroscopy 1988 Vol. 42 No. 5 p. 865. 

Continuum regression comprises an infinite-member family of methods for multivariate calibration. PLS 
and PCR are individual members of the continuum regression family. See, M. Stone, and R. J. Brooks 
(1990), "Continuum Regression: Cross-validated Sequentially Constructed Prediction Embracing Ordinary 
Least Squares, Partial Least Squares and Principal Components Regression," Journal of the Royal 

55 Statistical Society B., 52, pp. 237-269. 

Two or more analytes can be calibrated or analyzed simultaneously by using a global PLS method 
called PLS2. In practice PLS2 tends to underperform PLS1 (analytes calibrated sequentially) primarily 
because model complexity is fixed for all analytes. See, Lindberg W. Pesson, J. A. Wold, S. (1983), 
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Analytical Chemistry, 55, p. 643. 

Ridge regression is another multivariate calibration method which has been used in non-medical 
situations in which the intensities at different spectral frequencies exhibit significant collinearity and the 
number of calibration samples exceeds the number of spectral frequencies. Hoerl et al., Practical use of 

5 ridge regression: a challenge met, "Applied Statistics" 34, 114-120, (1985), showed that ridge regression 
was a viable competitor to multiple linear regression in the context of predicting percent protein in wheat 
samples by using reflectance in the near Infrared region. Naes et al.. Comparison of lineal statistical 
methods for calibration of NIR instruments, "Applied Statistics" 35, 195-206, (1985), concluded that ridge 
regression is a viable competitor to PLS and PCR when the number of spectral frequencies approaches the 

10 number of calibration samples. 

Another type of multivariate algorithm gaining wide acceptance is a pattern recognition technique using 
what are called neural networks. Like other multivariate calibration methods, neural networks learn from the 
input they are given. They have the potential advantage that they can explicitly model nonlinearitles. 
However, they also tend to be more susceptible to overfitting, and slower to compute, and are more difficult 

75 to interpret than PLS, PCR, and MLR. 

The applicants recognize that both the preferred embodiment of U.S. patent No, 4,975,581 and co- 
pending European patent application, Publication No. 0 522 674 A2 utilize the partial least square algorithm 
(PLS). However, the reasons for utilizing PLS and the other multivariate algorithms in the invention disclosed 
and claimed herein is quite different from the reasons it was utilized in the above described patent and 

20 application. For noninvasive glucose determination, the limiting factor in measurement is the lack of 
information available. The determination of a blood analyte, such as glucose, requires a very high signal-to- 
noise ratio and a sophisticated algorithm for extraction of a minuscule amount of information (glucose is, 
normally, 0.1 weight percent of blood). In the case of a pulse oximeter suitable for fetal monitoring, the 
information is abundant (i.e., babies that are profoundly hypoxic are blue), but the environment of operation 

26 is extreme. The reflected light-oxygen saturation relationship is highly nonlinear, the signal for analysis is 
extremely noisy and interfering background components must be removed by correlating with the pulsating 
blood. 

In contrast to the foregoing, the rationale for using multl-variate analysis in the present application of 
noninvasive blood gas determination is to enable accurate determination of blood gas parameters where the 

30 Information content in the spectral domain overlaps and where the infrared patterns for pH, PCO2. PO2 and 
[HCOs*] are small or do not exist in the absence of interactions with water and other blood or tissue 
components. Examination of Figure 10 shows that the regions of spectral information for the various 
components overlap extensively. This is especially true for pH which does not exhibit a strong correlation in 
any specific region. The problem of differentiating bicarbonate from carbon dioxide is also difficult as 

35 evidenced by their similar correlation curves. 

The use of visible and infrared spectroscopy for the noninvasive determination of blood gas parameters 
is not obvious. In fact "at first glance, the determination of pH by infrared spectroscopy is so implausible as 
to seen ridiculous." "Salinity determination using NIRH." T. Hirschfeld, Applied Spectroscopy , Volume 39, 
Number 4, 1985. Not only does water have exceptionally strong absorption bands, but H+ and O2 sat. have 

40 no absorption bands of their own in the near infrared. 

Infrared spectroscopy obtains quantitative and qualitative infonnation from the vibrational motion in 
molecules. Each type of chemical bond in a molecule will absorb different frequencies of infrared energy, 
giving rise to the characteristic patterns seen in an infrared spectrum. These patterns, called absorption 
bands, change in magnitude with concentration. Traditional methods of quantitative spectroscopy correlate 

45 band height or the area of a band to the concentration of the species under study. In the present application 
of infrared spectroscopy, characteristic infrared patterns for each species under study (pH, 02, CO2. 
[HCO3"] are small or do not exist when not interacting with other components. However, these species do 
have an effect on species that do absorb, specifically hemoglobin and water. Hydrogen ion, being an ion 
rather than a molecule, does not have infrared bands. However, hydrogen ions will bind to other species in 

50 solution that are infrared active, thus a correlation for pH can be based on secondary spectroscopic effects. 
Oxygen, although a molecular species, does not have Infrared bands. However, It will also interact with 
other species that are infrared active and produce a change in the spectra. In particular, oxygen reacts with 
hemoglobin, producing a marked change in the infrared spectra of hemoglobin between 600 and 1000 nm. 
This reaction is nonlinear, however, and quantitative models must account for this nonlinearity. For CO2. 

55 and [HCOa~], near infrared spectral absorbances are small. CO2 can be quantitated easily in the mid- 
infrared region. [HCOs"] also has bands In the mid infrared region. But In the near Infrared region, CO2 and 
[HCOa"] absorbance are small compared to the much larger effect of hemoglobin and proteins. Again, 
secondary effects must be studied. Thus, accurate measurement of the various blood gas parameters is 
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difficult due to the small absorbances and/or secondary effects resulting from concentration changes. These 
differences necessitate the use of multivariate analysis. 

The clinical utility of a noninvasive blood gas monitor, (for the simultaneous determination of pH, PCO2, 
PO2. [HCOa''] and O2 saturation) Is clearly desirable and the technology for realization of such a monitor as 
5 well as data utilizing the technology is disclosed herein. Continuous monitoring would be of significant 
benefit when a patient's cardiopulmonary status is changing rapidly. 

It is an object of the present invention to provide an instrument which will represent a remarkable 
improvement in blood gas monitoring technology by: 

A. Providing noninvasive blood gas monitoring which can determine all blood gas parameters in arterial 
10 blood and/or tissue. 

B. Providing continuous determination of ail blood gases in real time; and 

C. Having the ability to provide the physician with a measure of validity or assurance of accuracy by 
employing outlier detection methods. 

The ability to determine inaccurate results is extremely important, especially when caring for the 
15 critically ill patient. 

It is an object of the present invention to determine the blood gas parameters in either arterial blood or 
human tissue. In the case of tissue determination, the optical determination will measure the blood gases of 
all the components that are irradiated by the light. Although total tissue determination is not currently used It 
may be a better measure than arterial blood gas measurement because the physician is interested In 

20 knowing how well the tissue is being perfused. 

It is an additional object of the present invention to determine blood gas parameters in arterial blood. 
The method of sampling during the systolic and diastolic portions of the cardiac cycle enables determina- 
tion of the pulse blood spectra. This spectra can be subsequently analyzed by multivariate algorithms for 
determination of blood gas parameters. Thus, the present invention enables measurement of blood gas 

26 parameters in arterial blood noninvasively. 

In the case of the present invention, it is an object to use wavelengths in the 500-2400 nm range to 
determine arterial blood gas parameters. As several components are determined, measured intensities at 
several wavelengths will be required. Thus, determination of multiple components may require measure- 
ment of three or more wavelengths. 

30 An additional advantage of the methodology of the present invention is the ability to report all blood gas 
parameters. The various parameters can be measured directly or calculated from well established equa- 
tions. The relationship between pH, PCO2, and [HCO3"] is described by the Henderson-Hasselbach 
equation: 

35 

pH = 6.1 + log [HCOJ] 

0.03 * PCO2 

40 The Henderson-Hasselbach equation allows for determination of two components with calculation of the 
third. 

Additionally PO2 can be calculated from O2 sat. and vice versa. O'Riordan et al. compared various 
calculation algorithms and assessed their accuracy relative to normal human data. See. J. F. O'Riordan. T. 
K. Goldstick, L N. Vida, G. R. Honig, and J. T. Ernest, "Modelling whole blood oxygen equilibrium; 
45 comparison of nine different models fitted to normal human data", Advances in Experimental Medicine, Vol 
191. pp 505-522, 1985. The concurrent use of pH In these calculations to more accurately assess the Pso 
value may improve overall calculation results. 

Summary of the Invention 

50 

This invention relates to methods and apparatus for, preferably, determining non-invasively and in-vivo 
at least two of the five blood gas parameters (i.e., pH, PO2, PCO2, [HCO3"] and O2 sat.) in a human. 
The preferred method includes the steps of: 

a. generating two or more different wavelengths of light, wherein the wavelengths are in the range of 500 
65 nm to 2,400 nm; 

b. irradiating in vivo and non-invasively blood containing tissue with such wavelengths, so that there Is 
differential attenuation of at least some intensities of the wavelengths by the blood containing tissue as a 
function of the wavelengths, the differential attenuation causing intensity variations of those wavelengths 
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incident from the blood containing tissue as a function of those wavelengths, the tissue and the unknown 
levels of pH, [HCO3*"]. PO2 and O2 sat.; 

c. measuring said intensity variations from the blood containing tissue to obtain a set of intensity 
variations v. wavelengths: 

5 d. generating at least one calibration model which is a function of conventionally measured blood gas 
parameters obtained from a set of calibration samples and intensity variations v. wavelengths data 
obtained from irradiating such set of calibration samples with the same 2 or more different wavelengths 
of light; and 

e. calculating the values of the unknown blood gas parameters in the blood containing tissue from the 
70 measured Intensity variations from the tissue utilizing at least one algorithm and the calibration model, 
the calculation using 2 or more variables. 

The method and apparatus: (1) provide for the detection of both spectral and contractration outliers in 
both the blood containing tissue and the calibration set; (2) provide for the determination of the blood gas 
parameters in both tissue and arterial blood; (3) utilize one or more algorithms, at least one of which is a 
76 multivariate algorithm such as PLS, PLR and neural networks; (4) includes the pretreatment of data; (5) 
utilizes a calibration model for each parameter being detemnined; and (6) can be used with both reflected 
and transmitted light. 

The apparatus includes, for instance, a spectrometer apparatus for conducting light to and from the 
patient, an array detector, a computing unit, and a visual display module. 
20 The method and apparatus can also be used invasively. 

Brief Description of the Drawings 

Figure 1 is a schematic illustration of the preferred embodiment of the invention; 
25 Figure 2 is a graphical representation of the basic principle of how additional blood spectra can be 
determined; 

Figure 3 is a series of graphs comparing univariate calibration to multivariate calibration; 
Figure 4 is a chart showing the general approach used in multivariate statistical methods to generate a 
mathematical calibration model and to use this model to quantitatively determine concentrations and/or 
30 properties from the spectra of unknown samples; 

Figure 5 is a schematic of the test apparatus used in the sheep studies; 

Figure 6 is a graph of the electrical activity of the heart, the arterial pulse pressure, and the sampling 
period and corresponding electrical triggers; 

Figure 7 is a graph illustrating the distribution of [HCO3"] and pH determinations obtained with tests 
35 utilizing the apparatus illustrated in Figure 5; 

Figure 8 is a scatterplot matrix illustrating the relationships between reference [HC03~], O2 sat., PCO2, 

pH and PO2 obtained during the sheep studies; 

Figure 9 displays tissue spectra from several representative samples; 

Figure 10 Is a graph illustrating the relationship between the absolute value of the correlation coefficient 
40 and wavelength for each of the five blood gas parameters, 

Figure 11 is a plot of predicted tissue pH (obtained from the tests utilizing the test apparatus illustrated in 

Figure 7) vs. measured arterial blood pH using the best ratio algorithm; 

Figure 12 is a plot of predicated tissue pH vs. measured arterial blood pH utilizing MLR; 

Figure 13 is a plot of predicted tissue pH vs, measured arterial blood pH utilizing PLS; 
45 Figure 14 is a plot of predicted tissue pH vs. measured arterial blood pH utilizing the neural networks 

algorithm: 

Figure 15 is a plot of predicted tissue [HCOs^] vs. measured arterial blood [HCOs^] utilizing the best 

ratio algorithm; 

Figure 16 is a plot of predicted tissue [HCOa"] vs. measured arterial blood [HCOa"] utilizing MLR; 
50 Figure 17 is a plot of predicted tissue [HCOa"] vs. measured arterial blood [HCOa"] utilizing PLS; 

Figure 18 is a plot of predicted tissue [HCOa'] vs. measured arterial blood [HCOa"] utilizing the neural 

networks algorithm; 

Figure 19 is a plot of predicted tissue PCO2 vs. measured arterial blood PCO2 utilizing the best ratio 
algorithm; 

55 Figure 20 is a plot of predicted tissue PCO2 vs. measured arterial blood PC02 utilizing MLR; 
Figure 21 is a plot of predicted tissue PCO2 vs. measured arterial blood PCOa utilizing PLS; 
Figure 22 is a plot of predicted tissue PCO2 vs. measured arterial blood PCO2 utilizing the neural 
networks algorithm; 
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Figure 23 is a plot of predicted tissue PO2 vs. measured arterial blood PO2 utilizing the best ratio 
algorithm: 

Figure 24 is a plot of predicted tissue PO2 vs. measured arterial blood PO2 utilizing MLR; 
Figure 25 is a plot of predicted tissue PO2 vs. measured arterial blood PO2 utilizing PLS; 
5 Rgure 26 is a plot of predicted tissue PO2 vs. measured arterial blood PO2 utilizing the neural networks 

algorithm; 

Figure 27 is a plot of predicted tissue O2 sat. vs. measured arterial blood O2 sat, utilizing the best ratio 
algorithm; 

Figure 28 is a plot of predicted tissue O2 sat. vs. arterial blood measured O2 sat. utilizing MLR; 
10 Figure 29 is a plot of predicted tissue O2 sat. vs. arterial blood measured O2 sat. utilizing PLS; 

Figure 30 is a plot of predicted tissue O2 sat. vs. arterial blood measured O2 sat. utilizing the neural 
networks algorithm; 

Rgure 31 is a graph of several systolic and diastolic spectra obtained by EKG Lock; 
Figure 32 is a graph of the optical pulse pressure variation during EKG data acquisition; 
15 Rgure 33 is a plot of a couple of systolic and diastolic spectra obtained by EKG lock following 
normalization; 

Figure 34 is a graph of the optical pulse pressure variation of the normalized data taken by EKG lock 
acquisition: 

Figure 35 is a plot of a representative pulse blood spectrum; 
20 Figure 36 is a plot of predicted pulse blood spectrum pH (obtained from the tests utilizing the test 

apparatus illustrated in Figure 7) vs. measured arterial blood pH using the best ratio algorithm; 

Figure 37 is a plot of predicted pulse blood spectrum pH vs. measured arterial blood pH utilizing MLR; 

Rgure 38 is a plot of predicted pulse blood spectrum pH vs. measured arterial blood pH utilizing PLS; 

Figure 39 is a plot of predicted pulse blood spectrum [HCOs"] vs. measured arterial blood [HCO3'] 
25 utilizing the best ratio algorithm; 

Figure 40 is a plot of predicted pulse blood spectrum [HCOa*] vs. measured arterial blood [HCO3"] 

utilizing MLR; 

Figure 41 is a plot of predicted pulse blood spectrum [HCO3"] vs. measured arterial blood [HCO3"] 

utilizing PLS; 

30 Figure 42 is a plot of predicted pulse blood spectrum PO2 vs. measured arterial blood PO2 utilizing the 
best ratio algorithm; 

Figure 43 is a plot of predicted pulse blood spectrum PO2 vs. measured arterial blood PO2 utilizing MLR; 
Rgure 44 is a plot of predicted pulse blood spectrum PO2 vs. measured arterial blood PQ2 utilizing PLS; 
Figure 45 is a plot of predicted pulse blood spectrum O2 sat. vs. arterial blood measured O2 sat. utilizing 
35 the best ratio algorithm; 

Figure 46 is a plot of predicted pulse blood spectrum O2 sat. vs. arterial blood measured O2 sat. utilizing 
MLR: 

Figure 47 is a plot of predicted pulse blood spectrum O2 sat. vs. arterial blood measured O2 sat. utilizing 
PLS: 

40 Figure 48 is an alternate schematic of a noninvasive ABG monitor using multiple sources and a single 
detector: 

Figure 49 is another alternate schematic of a noninvasive ABG monitor using a single source and 
multiple detectors; and 

Figure 50 is a schematic of an invasive ABG monitor of the present invention. 

45 

Description of the Preferred Embodiment 

With reference to Figure 1, the preferred noninvasive blood gas monitor 11 of the present invention 
includes a spectrometer 13, an electronics and computer processing module 15, and a visual display 

50 module 17. Spectrometer 13 includes a broad band tungsten halogen light source 19, a focusing mirror 21, 
a fiber optic interface 23, a second fiber optic interface 25, a grating 27. an array detector 29. and an 
electronic buss 31. Module 15 includes a computing unit which, in turn, contains a microprocessor, 
memory, the data preprocessing algorithms, the multivariate calibration model, the multivariate algorithm 
and outlier detection algorithms. Visual display module 17 includes a pH display 41, a partial pressure of 

55 carbon dioxide display 43, a bicarbonate display 45, a partial pressure of oxygen display 47, and an oxygen 
saturation display 49. As illustrated, each of these displays includes each determination and the estimated 
uncertainty of such determination. Also Included in module 17 are a pH trend display 51, a partial pressure 
of carbon dioxide trend display 53, a bicarbonate trend display 57, a partial pressure of oxygen trend 
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display 57, and a oxygen saturation trend display 59. Visual display 17 also includes a tissue/pulse mode 
switch 60, which permits the doctor or other operator to change the mode of instrument operation from the 
tissue mode to the blood pulse mode, or vice-versa. Finally, display module 17 includes lights 61, 63, 65, 
67 and 69 to indicate if the associated blood gas parameter which has just been determined is an outlier 
5 (e.g.. light 61 is energized when the pH determination is an outlier). 

To transmit light from spectrometer 13 to the fingertip 71 of the patient being monitored, monitor 11 
includes a source fiber optic bundle 73, which terminates at finger/fiber device 75. Receiving fiber optic 77 
returns the light from finger/fiber holder 75 to fiber optic interface 25. Finger/fiber holder 75 allows 
transmission through finger 71 and allows for adequate holding of the finger. 

70 In operation, source 19 emits selected frequencies from approximately 500 to 2400 nm. This light is 
focused on the end of fiber optic 73 held in interface 23 via focusing mirror 21 and then transmitted via 
source fiber 73 to illuminate the tissue, bone, nail and blood in fingertip 71 , The portion of the light which is 
transmitted through fingertip 71 is then returned to spectrometer 13 by fiber bundle 77. The returning light 
is then separated into various wavelengths and detected by the array detector 29, capable of detecting at 

75 least some of the wavelengths of light between 500 to 2400 nm. 

The reflected light intensities at the various frequencies are then analyzed by the microprocessor of 
module 15 employing, a multivariate algorithm (such as PLS or PGR) utilizing several wavelengths from 
various regions of the entire spectral range of the transmitted light. Instrument 11 can be operated in two 
different modes: (1) the tissue determination mode; or (2) the pulse blood mode. In the tissue determination 

20 mode of operation the spectral data from array detector 29 is preprocessed by appropriate algorithms in 
module 15. Such preprocessing to include trimming, averaging, scaling, normalization and derivatization as 
necessary. The resulting processed spectrum is then analyzed by the algorithm contained in module 15. In 
the prefen-ed embodiment the algorithm employed is partial least squares (PLS). Such an algorithm can 
either quantify directly from the spectra! data or through use of other equations calculate ail blood gas 

25 parameters. Only two of the following (pH. PC02 and [HCO's] need to be optically measured; with 
calculation of the third by the Henderson-Hasselbach equation. 

Instrument 11 will display the current values of pH, PCO2. [HCOa"], PO2 and O2 sat., their estimated 
uncertainties, as well as the past history of each parameter in displays 51. 53, 55. 57 and 59. The spectral 
data analyzed and the current determinations will be examined by outlier detection algorithms contained in 

30 the computing unit in module 15. If the spectra from the patient or the resulting determinations are 
abnormal then the instrument will indicate that the accuracy of those particular determinations are poor 
through the corresponding outlier lights, 61-69. A condition which might result in a spurious result might 
occur if the patient's fingertip were only partially inserted in the finger/fiber device 75. Other possibilities 
include the instrument being out of calibration, extreme blood gas values not included in the calibration 

35 model, or new and unique blood chemistries. 

In blood pulse mode, spectra are acquired throughout the cardiac cycle of the patient. The resulting 
spectral data are analyzed to establish which spectra correspond with maximum concentration of blood (or 
maximum dilation) In the arterial system of the patient, and which spectra correspond with minimum 
concentration (or minimum dilation) of the arterial system. The spectra associated with minimum dilation will 

40 contain information regarding blood, skin, bone, and tissue. The spectra associated with maximum dilation 
will contain the same information plus an additional amount of blood information. To determine optical pulse 
pressure, the spectral difference at the isobestic frequency is calculated. Absorbance by wavelengths in the 
region from 800 to 810 nm are not strongly influenced by changes in oxygen saturation and are proportional 
to the amount of hemoglobin present. Thus, a plot of the average value between 800-810 nm will show the 

45 variation in hemoglobin concentration. The magnitude of this difference will correspond to the spectral 
difference between the systolic and diastolic pulses. This spectral difference is, in turn, a general measure 
of the quality of the spectral data. Those spectra having an abnormal optical pulse pressure are trimmed 
(i.e., removed) prior to further analysis. The trim process removes spectra which deviate significantly from 
the norm due to motion of the patient or other artifacts. As some of these artifacts will influence the spectra 

60 in an adverse manner, the removal of these spectra is important for production of an accurate pulse 
monitor. Additional preprocessing of the data may be required due to specific operating conditions. The 
calculation of the "additional" blood spectrum can be performed by several methodologies, including 
subtraction of the appropriately transformed spectral data from the maximum and minimum dilation. The 
above process effectively removes the interfering background and provides the multivariate algorithm with a 

55 spectrum corresponding to the additional blood. The resulting additional blood or pulse blood spectra may 
be subsequently averaged with existing scans to improve the signal-to-noise ratio of the spectrum to be 
analyzed. The resulting spectrum is then analyzed by the multivariate algorithm contained module 15. In the 
preferred embodiment the algorithm employed would be partial least squares or neural networks. The 
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algorithm will either measure directly or calculate all blood gas parameters. The Instrument will display the 
current value of pH, PCO2, [HCOs"], PO2 and O2 sat., their uncertainties as well as the past history of each 
such parameter as explained above with regard to the tissue mode. The accuracy of each of these 
determinations will also be assessed through the use of outlier detection methods contained in the 

5 computing unit in module 15. 

An additional embodiment of the invention includes apparatus (not shown) for obtaining information 
regarding the electrical activity of the patient's heart, which can provide information to assist in determina- 
tion of maximal and minimal arterial dilation. Typically most patients in intensive care units will have an EKG 
monitor in operation and use of this early accessible information will assist in the operation of Instrument 1 1 

10 In the pulse mode. With reference to Figure 6, maximum expansion of the arterial system due to ventricular 
contraction occurs at a set interval following the R peak of the QRS complex. The QRS complex is created 
by the depolarization of the ventricular muscle. This complex precedes ventricular contraction which results 
in ejection of blood from the heart. Minimum expansion of the arterial system Is present prior to ventricular 
contraction and corresponds to a time period in the vicinity of the P-wave. The P-wave results from 

75 depolarization of the atrial muscle. The time of minimum expansion is related to both the R peak of the 
QRS complex and heart rate. Correlation with the electrical activity of the heart may be desired for 
simplification in the processing of the spectral data and may be necessary for effective operation during 
period of decreased pulse pressure or excessive finger movement. To summarize, the electrical activity of 
the patient's heart provides additional information for Improved operation, especially under adverse 

20 conditions. 

It is the authors* experience that pretreatment of the spectral or concentration data often times improves 
the analysis precision in the calibration set and unknown analyses, as well as increases the robustness of 
the multivariate calibration models. Thus, data pretreatments including but not limited to centering, scaling, 
normalizing, taking first or higher order derivatives, smoothing. Fourier transforming, and/or linearization can 

26 all improve the analysis precision and accuracy. Additionally, trimming of those spectra differing from the 
norm for any reason can improve overall instrument performance. Such trimming techniques are not limited 
to assessment of the optical pulse. Other techniques include the comparison of the summation of the area 
under the curve, comparison of individual spectra with the average spectra, comparison of the Fourier 
transform, examination of the characteristic modes of vibration in the spectra, and comparison of principle 

30 component regression factors. These various techniques may be used individually or together for effective 
trimming. These pretreatments can also improve the robustness of the model to instrument drift and can 
improve the transfer of the calibration model between instruments. 

To demonstrate the technical feasibility of the invention, noninvasive determination of both tissue and 
arterial blood gas parameters were performed on baby lambs under a variety of physiological conditions. 

35 The object of the investigation was to induce all reasonable acid-base disturbances in the lambs. 
Specifically, the acid-base balance of the lambs was manipulated so as to induce respiratory acidosis, 
respiratory alkalosis, metabolic acidosis and metabolic alkalosis. 

Newborn lambs (age 10-25 days) were instrumented as schematically illustrated in Figure 5. The 
animal's respiration was controlled by a conventional mechanical ventilator 81 with endotracheal intubation 

40 83. Metabolic physiology was controlled by administration of bicarbonate or lactic acid as appropriate 
through a catheter 85. Cardiac physiology was monitored by placement of arterial and venous trans- 
ducer/monitor and external placement of standard EKG electrodes 89 connected to EKG monitor 91. 
Oxygen saturation was monitored by use of a standard pulse oximeter 93. 

After Intubation and instrumentation of the lamb as set forth above, a spectrometer 101 equipped with a 

45 source fiber 103 and a receiving fiber 105 was applied to the lamb for procurement of spectral data. The 
exit end of source fiber 103 was positioned on one side of and in contact with the lamb's shaved leg, with 
receiving fiber 105 positioned on and in contact with the opposite side of the leg. Light from a tungsten 
halogen light source 107, covering the frequency range 500-1000 nm, was coupled to fiber 103 and 
transmitted to the lamb's leg. The light exits fiber 103 and enters the lamb's tissue where it interacts with all 

50 components of the lamb's (e.g., blood, interstitial fluid, muscle, bone, fat and skin), before exiting from the 
opposite side of the leg. At this point some of this light enters the receiving fiber optic. The light entering 
into receiving fiber 105 is returned to spectrometer 101 where it is separated by wavelength by using a 
holographic grating and subsequently detected on a silicon array detector 109. The light intensities at 
detector 109 is then converted to digital numbers and stored in computer 117. 

55 During the lamb study, spectra, EKG information, and the arterial pulse were all acquired simulta- 
neously. Computer 115 is used to process EKG signals from monitor 91 and subsequently to trigger 
spectrometer 101, via trigger line 111. Computer 117 is used to record spectral data from spectrometer 101 
via electronic buss 113. The spectral data were acquired in two different modes: (1) pulse gated; and (2) 
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EKG lock. In the pulse gated mode the spectral data was acquired at the same point in the cardiac cycle. 
Specifically, the shutter (not shown) to silicon array detector 109 was opened with detection of the R wave 
of the QRS complex of the EKG signal. Detector 109 was exposed to the light transmitted through the 
lambs leg for 50 msec, after which the shutter was closed. The spectral data was taken in this manner so 

5 the same or approximately the same amount of blood was present in the tissue during the spectral 
acquisition period. The data sampled in this manner is subsequently referred to as "pulse gated data". The 
light detected in pulse gated mode has interacted with the tissue as explained above and contains 
information from the cube of tissue through which it passed. The pulse gated data was then subsequently 
analyzed by several algorithms to predict arterial blood gas parameters. 

10 The second mode of spectral acquisition is referred to as EKG lock. The concept is to use the electrical 
activity of the heart to determine the correct time to obtain systolic and diastolic spectra. The time 
separation between two consecutive R waves of the QRS complex was determined and the heart rate 
subsequently calculated. Spectral acquisition corresponding to the diastolic portion of the cardiac cycle 
occurred at a fixed time interval of 30 msec after the R wave. The systolic spectral sampling occurred at 

75 50% of the cardiac cycle. The resulting spectral data is stored on computer disk as previously described. 

Generation of Non-Physiological Acid-Base Conditions 

The acid-base status of the lamb was varied over the maximum range reasonable. To introduce 
20 respiratory perturbations the partial pressure of carbon dioxide in the lamb's blood was manipulated by 
changing the rate of mechanical ventilator, particularly the minute ventilation to the lamb. To introduce 
metabolic perturbations the lamb was given intravenous bicarbonate or lactic acid. Additionally, to examine 
the ability of the Instrumentation to determine both O2 sat. and PO2, the fraction of inspired oxygen (F1O2) 
was varied randomly over the course of the tests. The animal's O2 sat. was varied from 80 to 100% while 
25 the partial pressure of oxygen varied from 45 to 115 mm Hg. O2 sat. and PO2 are related quantities, but are 
not related in a linear fashion. 

The arterial blood samples from the lambs were measured on a Mailinckrodt Sensor System Gem-Stat 
blood gas/electrolyte analyzer. The Gem-Stat is a commonly used bedside blood gas monitor and is a 
clinically accepted method for blood gas determination. 
30 The accuracy and precision of the reference values used in this study were assessed by analyzing 
summary data provided by Zaioga et al. (1989). Based on Zaioga et al., our estimates of the analyte- 
specific standard deviation are as follows: 
pH - .025 

PO2 - 10 to 20 mm (10 at 50mm PO2 and 20 at 120 nm PO2) 
35 PCO2 - 3mm 

[HCOa"] - 2 millimoles per liter 

The accuracy of the measured O2 saturation, performed on a Radiometer, 0SM2, was difficult to 
assess due to lack of information available from the manufacturer. Nevertheless, the authors' experience 
with the Hemoximeter suggests that the actual standard deviation of repeated measurements is approxi- 
40 mately 2%. 

Figure 7 shows reasonable physiological k)oundaries for pH. bicarbonate, and carbon dioxide. One of 
the objects of the lamb study was to map out the entire physiological region, of respiratory acidosis and 
alkalosis as well as metabolic acidosis and alkalosis. Examination of Figure 7 shows the acid-base balance 
of the lamb was disturbed over the entire physiologically reasonable space. Scatter plot, Figure 8, shows 

45 the variation in oxygen saturation and partial pressure of oxygen over the course of the study. As pH, PCO2 
and [HC03~] are all constrained by the Henderson-Hasselbach equation, it is essential to make sure that 
the variation in one variable does not correlate completely with the variation in another variable. 

Figure 8 is a scatterplot matrix showing the relationship between all blood gas parameters. Examination 
of the figure shows that PO2 and O2 sat. are correlated as one would expect. However, there does not 

50 appear to be significant correlation between pH, PGO2 and [HCOa']. This lack of correlation demonstrates 
that each parameter (other than O2 sat. and PO2) was varying indeF)endently. 

Analysis of Spectral Data in Simulated Tissue Mode. 

55 The pulse gated data was acquired as stated above through the lamb's leg. Thirty scans were taken 
and subsequently averaged to form a single spectrum for each sample or observation point. Several 
representative spectra are plotted in Figure 9. These raw intensity data were subsequently processed and 
analyzed in several ways. 
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As previously stated the optica! measurement of all arterial blood gas parameters is difficult, as the 
spectral information for a given analyte typically overlaps with information regions for other analytes. The 
overlapping of spectral information makes it impossible to make accurate determinations of blood gas 
parameters by univariate algorithms, if the concentration of the overlapping analyte changes and the 

5 magnitude of the spectral changes Introduced by the overlapping analyte Is similar to that of the analyte of 
interest. Figure 10 shows the relationship between the absolute value of the correlation coefficient and 
wavelength for each blood gas parameter. Examination of Figure 10 clearly shows that the regions of 
spectral information overlap. 

To demonstrate that various algorithms can be used to measure some or all blood gas parameters, the 

70 following algorithms were investigated: (1) best ratio analysis; (2) multiple linear regression: (3) partial least 
squares: and (4) neural networks. 

In comparing the capabilities of the various algorithms it is important to note that the standard error of 
calibration (SEC) and standard error of prediction (SEP) are presented depending upon the analyte. sample 
spectra, or algorithm used. These two measures of predictive ability are not the same and have subtle but 

75 important differences. The SEC refers to errors in predicting sample analyte values that were used to build 
the model. The SEP is reported if a cross validation technique was used. When all samples are 
simultaneously used to develop the calibration model, the SEC reported can be less than the standard error 
of the reference. Such an occurrence suggests overfitting of the experimental data by the calibration model. 
The more sophisticated algorithms (e.g. PLS, PCR and neural networks) are all susceptible to some 

20 overfitting of this nature when all calibration samples are simultaneously included in the calibration model. If 
a cross validation technique is used, the standard error of prediction (SEP) can be reported. In the cross 
validation procedure, one or more samples are omitted from the calibration, and the calibration model is 
determined based upon this reduced sample set. The concentration of the omitted sample(s) is then 
predicted using the predetermined calibration model. The process is repeated until each sample has been 

26 left out of the calibration model. The SEP, based on cross validation provides a realistic measure of the 
instrument's prediction ability in the clinical/patient setting. In the discussion that follows, SEC is reported in 
some cases, while SEP is reported in others. The errors are reported differently due to available computer 
software, and computational time required. In addition, the errors present in the reference data cleariy limit 
the potential precision of any model. Many of the figures of predicted vs. reference concentrations show 

30 large apparent scatter. Often times, this scatter is a direct result of the limited precision of the reference 
methods. When the prediction precision is the same as the precision of the reference method, then the 
plots will have the same appearance as if a set of samples were measured two times by the reference 
method and the reference value of the first measurement were plotted vs. the second. 

35 Multiple Linear Regression 

Multiple Linear Regression was used to analyze the Gaussian averaged intensity spectra. The correla- 
tion coefficient between the spectral intensity at each wavelength and the analyte concentration were 
calculated. The wavelength exhibiting the best correlation coefficient was chosen from the 384 data points 

40 covering the wavelength range 640-970 nm as the first wavelength in the model. Subsequent wavelengths 
were chosen based on correlations between the residuals which remained after information from the 
previously chosen wavelength(s) was removed. The optimal number of intensities for analysis was based on 
the model that yielded the minimum in the cross-calibration SEP. 

Table 2 lists the parameter, number of frequencies, the number of samples analyzed, and the standard 

45 errors of prediction. The results of the analysis are illustrated in Figures 12, 16, 20. and 24. The standard 
errors of prediction were calculated using a cross validation procedure leaving one sample out at a time. 
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TABLE 2 



10 



MLR - Lamb Tissue 


Parameter 


Number of 
Wavelengths used 


Number of Samples 


SEP 


PH 


2 


39 


0.07 


PCO2 


2 


38 


4.9 


PO2 


1 


36 


6.9 


[HCO3-] 


3 


40 


2.7 


O2 sat. 


1 


40 


2.6 



75 

Partial Least Squares 

Partial Least Squares was used to analyze the second derivative of the raw intensity spectra. A 10 point 
20 moving average followed by two first difference calculations was used to calculate the derivative. A separate 
PLS model was formed for each of the five parameters. Table 3 lists the parameters, the number of PLS 
factors, the wavelength range, the number of samples, and the standard errors of prediction. The standard 
errors of prediction were calculated using a cross validation procedure. The results of the analysis are 
illustrated in Figures 13. 17, 21, 25, and 29. 

25 

TABLE 3 



PLS - l-amb Tissue 


Parameter 


Number of 
Factors 


Wavelength Range 


Number of Samples 


SEP 


pH 


7 


700-800nm 


39 


0.03 


PCO2 


1 


652-955nm 


38 


4,4 


PO2 


3 


652-955nm 


36 


6.0 


[HCO3-] 


3 


652-955nm 


40 


2.6 


O2 sat. 


3 


640-970nm 


40 


2.7 



40 

Neural Networks 

A neural network was developed for each of the parameters using a subset of the calibration samples. 

45 The subset was chosen randomly Table 4 lists the parameter, the number of samples analyzed, and the 
standard errors of the calibration. For all analytes measured, the reported SEC Is better than the precision 
of the reference method. This condition indicates that the neural network model relating spectral and 
concentration information has overfitted the data, and the results are unrealistically good. Nevertheless, 
neural network analysis of the data worked exceptionally well and the algorithms performance indicates that 

50 it is a functional algorithm for the present invention. Neural network cross validation was not performed 
since the cross-validation algorithms were not available and the computation time required is enormous. 
The results of the analysis are illustrated in Rgures 14, 18, 22, 26, and 30. 
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Neural Network - Lamb Tissue 


Parameter 


Number of Samples 


SEC 


pH 


34 


0.01 


PCO2 


31 


2.3 


PO2 


29 


2.1 


IHCO3-] 


35 


1.4 


02 sat. 


35 


0.4 



76 

Analysis of Spectral Data fn Simulated Puise Blood Mode 

The EKG lock data was acquired as stated above tiirough the lamb's leg. Twenty five systolic and 
twenty five diastolic scans were taken at each blood sample observation. Figure 31 shows five systolic and 

20 five diastolic spectra from the lamb's leg. Careful examination of the scans shows some difference between 
the systolic and diastolic scans. To quantify the difference between the systolic and diastolic spectra, the 
optical pulse pressure was calculated. The optical pulse pressure was determined by averaging the 
intensity values from 800 to 810 nm for each spectra. This region of the spectrum was chosen as it Is an 
isobestic region with respect to oxygen saturation. Thus, intensity measurements in this region of the 

25 spectrum vary with changing hemoglobin concentration and are not influenced by changes in oxygen 
saturation. As changing hemoglobin concentration will be proportional to changes in arterial expansion and 
contraction, analysis of the Intensity fluctuations In this region enables quantification of the optical pulse 
pressure. Rgure 32 is a plot of the optical pulse calculated from Intensity values at 800 to 810 nm as a 
function of scan number. 

30 Examination of Figure 32 reveals several important characteristics: (1) the optical pulse pressure is 
quite small; and (2) the summed isobestic intensity values during both systole or diastole are not constant 
during the sampling period. Variations in the intensity of light transmitted (i.e. optical puise pressure) arise 
from changes in the blood fractional volume of the tissue which alters light absorption. Changes in blood 
fractional volume are correlated with the capillary density in the dermis. The end result is that due to low 

35 capillary densities the optical pulse is very small. The small optical pulse makes pulse blood measurements 
from the leg of the lamb extremely difficult due the minuscule difference between the systolic and diastolic 
spectra. 

The problem of procuring a pulse blood spectrum is compounded by the variations In the spectra as 
shown in Figure 31 and in the variations observed in either the systolic or diastolic Isobestic sum. as shown 

40 in Figure 32. Examination of all fifty spectra indicated that some spectra were clearly abnormal while some 
simply had a baseline difference. To procure the best possible pulse blood spectra, the fifty samples were 
trimmed and then normalized. 

The trimming process was based upon the optical pulse pressure observed between sequential systolic 
and diastolic pulses. The optical pulse pressure valves were ranked in lowest to highest order. Those 

45 systolic-diastolic scans which had a optical pulse pressure difference between the 40th and 90th percentile 
were considered to be representative scans. The samples selected are indicated by a filled square in Rgure 
32. Examination of these scans still revealed significant baseline differences. Thus, the selected scans were 
normalized by the sum of the intensity values over the entire spectrum. Following normalization, all systolic 
scans had the same intensity sum value over the entire spectrum. The diastolic scans were normalized to a 

50 different value, but the process was executed in the same manner. The normalization process can remove 
linear baselines without altering the shape of the spectrum. The same systolic and diastolic spectra shown 
in Figure 31 are, following normalization, plotted in Figure 33. Examination of Figure 33 reveals the majority 
of the initial differences between the systolic or diastolic groups was baseline in nature. Figure 34 is a plot 
of the resulting optical pulse pressure vs. scan number, following normalization of all systolic and diastolic 

55 scans. 

The normalized systoiic-dlastollc groups selected by the trim procedure were used to calculate the 
puise blood spectrum, as follows: 
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Pulse Blood Spectruin = log (diastolic) - log (systolic) 

( spectrum) ( spectrxm) 

The resulting scans were subsequently averaged to generate the pulse blood spectrum for each observa- 
tion. Figure 31 shows a plot of a pulse blood spectra obtained from the lamb's leg. The resulting spectra 
were subsequently analyzed by several different algorithms: best ratio; multiple linear regression; and partial 
least squares. As shown in the case of the tissue data, neural networks represent a powerful algorithm that 
is well suited to this type of analysis. However, analysis by neural networks is not included due to logistical 
problems in getting the data processed. 

The analysis of the pulse blood spectra was done in the same manner as the tissue spectra. The 
results of the analysis are shown in the following tables. The analysis of the spectra to quantify the partial 
pressure of carbon dioxide did not yield predictive relevance. Thus, no PCO2 predictions are Included. 

Best Ratio Analysis 



Best ratio analysis of the trimmed, normalized pulse blood spectra was done in the same manner as set 
forth above, including processing of the spectra with a moving Gaussian weighted average. The results are 
20 shown in Figures 36, 39, 42 and 45. 

TABLE 5 



25 



30 



Best Ratio - Pulse Blood 


Parameter 


Number of Samples 


SEC 


PH 


32 


0.05 


PO2 


42 


12.7 


[HCO3-] 


38 


3.65 


O2 sat. 


38 


3.55 



35 



40 



Multiple Linear Regression 

The results of MLR analysis of the pulse blood spectra is summarized in Table 6. The table lists the 
parameter number of wavelengths used in analysis and errors of calibration. The results are presented in 
Figures 37, 40, 43 and 46. 

TABLE 6 



45 



50 



MLR - Pulse Blood 


Parameter 


Number of 
Wavelengths used 


Number of Samples 


SEC 


pH 


7 


32 


0.03 


PO2 


1 


42 


9.6 


[HCO3-I 


4 


38 


2.4 


O2 sat. 


5 


38 


1.3 



55 
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Partial Least Squares 

The results of PLS analysis of the pulse blood spectra is summarized in Table 7. The table lists the 
parameter, number of factors, wavelength region, number of frequencies and standard error of calibration. 
5 The results are presented in Figures 38. 41. 44 and 47. 

TABLE 7 



70 



75 



PLS - Pulse Blood 


Parameter 


Number of 
Factors 


Wavelength Range 


Number of Samples 


SEC 


pH 


7 


700-800nm 


39 


0.02 


PO2 


3 


652-955nm 


36 


6.2 


[HCO3-I 


3 


652-955nm 


40 


0.7 


O2 sat. 


3 


640-970 nm 


40 


1.2 



20 



25 



30 



40 



45 



50 



55 



Although the data used to demonstrate proof of concept in the lamb study was recorded from 500 to 
1000 nm, this is not the only frequency region of interest. Specifically the frequency region from 1000 to 
2400 nm contains information on both hydrogen ion concentration and CO2. (E. Watson and E.H. 
Baughman, "On-line analysis of caustic streams by near-infrared spectroscopy", Spectroscopy, vol 2, no 1 
pp 44-48. The region from 1000 to 1400 nm is especially important as it allows for transmission 
measurements to be made easily through the finger, yet is not influenced by melanin content. Melanin 
which varies with ethnicity has a relatively flat absorbance response from 500 to 1000 nm, but ceases to 
absorb significantly above 1000 nm. 

The detector used in our sheep study was a silicon detector with a response range of approximately 
500 to 1000 nm. Detectors for measurement of frequencies in the region of 1000 to 2500 nm are numerous 
and include. Indium Gallium Arsenide, Indium Antimonide. Germanium, Lead Selenide, and Lead Sulfide. 
Additionally the range of response for silicon detectors has been increased to approximately 1200 nm 
through the use of coating powders and various techniques. 

With reference to Figure 48, a second embodiment of a noninvasive blood gas monitor 101 of the 
present invention includes a spectrometer 103, an electronics and computer processing module 105 and a 
visual display module 107. Spectrometer 103 includes several light emitting diodes 111-125, a focusing 
lenses 127, fiber optic interface 129, a second fiber optic interface 131, a focusing mirror 133, a detector 
135, and an electronic buss 137. l\^odule 105 includes a microprocessor, memory containing the data 
preprocessing algorithms, multivariate calibration model and algorithms, and outlier detection methods. 
Visual display module 107 includes a plH display 141. a partial pressure of carbon dioxide display 143, a 
bicarbonate display 145, a partial pressure of oxygen display, 147, and an oxygen saturation display 149. 
Also included in module 107 are a pH trend display 151, a partial pressure of carbon dioxide trend display, 
153, a bicarbonate trend display, 155, a partial pressure of oxygen trend display 157 and a oxygen 
saturation trend display 159. Visual display 107 also includes outlier lights 161-169 for each of the previous 
display analytes. Mode switch 171, which permits the doctor or other operator to change the mode of 
instrument operation from tissue mode to blood pulse mode, or vice-versa. 

To transmit light from spectrometer 103 to the fingertip 67 of the patient, whose blood gases are being 
monitored, monitor 101 Includes a source fiber optic bundle 173, which terminates at finger/fiber device 
175. Receiving fiber optic 177 returns the light from finger/fiber holder 175 to fiber optic housing 131. 
Finger/fiber holder 175 allows transmission through finger 179 and allows for adequate holding of the finger. 

In operation, multiple light emitting diode sources 111-125 emit light at discrete time intervals in the 
frequency region between 500 and 2400 nm. This light is focused on the end of fiber optic interface 129, 
via focusing lenses 127 and then via source fiber 173 to illuminate the tissue, bone, nail and blood therein. 
A portion of the light which is transmitted through fingertip 179 is then returned to spectrometer 103 by 
fiber bundle 177. The returning light is then focused by focusing mirror 133 onto detector 135, capable of 
detecting frequencies of light between 500 to 2400 nm. 

The light intensities from the various diodes are then analyzed by processing unit 105 employing a 
appropriate algorithm. The instrument 101 can be operated in two different modes controlled by mode 
switch 171: (1) tissue determination mode; or (2) pulse blood mode. 
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Instrument 101 will display the current values of pH, PCO2, [HCOa'], PO2 and O2 sat. as well as the 
past history of said analyte in said displays 151-159. If a given analyte determination indicates that the 
analysis might be unreliable the outlier lights 161-169 will indicate such to the operator. 

Present data Indicates that pH, PCO2, [HCO3'], PO2 and O2 sat. can be determined noninvasively 

5 using wavelengths between 500-2500 nm. In addition, information is expected in the longer wavelength 
range 2500-7500 nm. This region cannot be used for noninvasive transmission measurements due to the 
targe absorbances of water. However, use of wavelengths in this region are well suited for invasive, in-vivo 
measurements of blood gas parameters. Prior work by the authors has shown that glucose, urea and 
creatinine can be measured in biological systems by using attenuated total reflective sampling. This 

10 spectral region exhibits bands from the vibrational motion of molecules. [HCOa^] and CO2 have vibrational 
bands that can be used for calibration. Thus, use of mid-Infrared spectroscopy can be used to measure two 
blood gas parameters. 

Oxygen saturation and PO2 could be determined invastvely and in-vivo through the use of visible and 
near infrared spectroscopy. Specifically light transmitted by the fiber could be diffusely reflected by the 
16 blood and returned to the Instrument for subsequent analysis. Thus, the combination of mid infrared 
spectroscopy with its sensitivity for vibrational spectroscopy coupled with vislble/near-infrared spectroscopy 
wttii its sensitivity for electronic transition spectroscopy is ideally suited for invasive in-vivo continuous 
blood gas measurements. 

With reference to Figure 49, a third embodiment of a noninvasive blood gas monitor 201 of the present 
20 invention includes a spectrometer, 203, an electronics and computer processing module, 205, and a visual 
display module 207. Spectrometer 203 Includes a single broad band source 211, a focusing mirror 213, 
tinger holding device 215. diffuser 217, optical tilters 219-233, and detector bank 235 and an electronic 
buss, 237. Module 205 includes a microprocessor, memory containing the data preprocessing algorithms, 
multivariate calibration model and algorithms, and outiier detection methods. Visual display module 207 
25 includes a pH display 241, a partial pressure of carbon dioxide display 243, a bicarbonate display 245, a 
partial pressure of oxygen display, 247. and an oxygen saturation display 249. Also Included in module 207 
are a pH trend display 251. a partial pressure of carbon dioxide trend display, 253. a bicarbonate trend 
display, 255, a partial pressure of oxygen trend display 257 and a oxygen saturation trend display 259. 
Visual display 207 also includes outiier lights 261-269 for each of the previous display analytes. Mode 
30 switch 271, which permits the doctor or other operator to change the mode of instrument operation from 
tissue mode to blood pulse mode, or vice-versa. 

In operation, tight from source 21 1 is transmitted through tinger 281 held in holder 215. The transmitted 
light is diffused by diffuser 217, tittered by optical filters 219-233 and subsequentiy recorded by detector 
bank 235. 

35 The light intensities from the various detectors are then analyzed by processing unit 205 employing a 
appropriate algorithm. The instrument 201 can be operated in two different modes controlled by mode 
switch 271: (1) tissue determination mode; or (2) pulse blood mode. 

Instrument 210 wilt display the current values of pH, PC02. [HCOa"], PO2 and O2 sat. as well as the 
past history of said analyte in said displays 251-259. If a given analyte determination indicates that the 

40 analysis might be unreliable the outiier lights 261-269 wilt indicate such to the operator. 

With reference to Figure 50, invasive blood gas monitoring can be performed by the Instrument shown 
In 301, the instrument includes a spectrometer, 303, an electronics and computer processing module, 305, 
and a visual display module 307. Spectrometer 303 includes a broad band light source 311 a focusing 
mirror 313, fiber optic interface 315, a second fiber optic interface 317, a grating 319, an array detector 321, 

45 and an electronic buss, 323. Module 305 includes a microprocessor, memory 35 containing the data 
preprocessing algorithms, multivariate calibration model and algorithms, and outiier detection methods. 
Visual display module 307 includes a pH display 341. a partial pressure of carbon dioxide display 343, a 
bicarbonate display 345, a partial pressure of oxygen display, 347, and an oxygen saturation display 349. 
Also included in module 307 are a pH trend display 351 , a partial pressure of carbon dioxide trend display, 

50 353, a bicarbonate trend display, 355, a partial pressure of oxygen trend display 357 and a oxygen 
saturation trend display 359. Visual display 307 also includes outlier lights 361-369 for each of the previous 
display analytes. 

To transmit light from spectrometer 303 to the arterial vessel 361, in which blood gases are being 
monitored, monitor 301 includes a source fiber optic bundle 373, which terminates at junction 375 with 
55 subsequent formation of a source and receiving fiber 377. The receiving portion of optical bundle 377 
returns the light to junction 375 with subsequent return of light to the instrument occurring through receiving 
fiber optic 379. 
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In operation, light source 311 or sources (not shown) emit light in the visible or infrared frequency 
regions. This light is focused on the end of fiber optic interface 315, via focusing mirror 313 and then via 
fibers 373 and 377 the light illuminates blood at the end of the fiber 363. The fiber optic probe is inserted 
through the skin 381 and into arterial vessel 361 . A portion of the light which has interacted with the blood 

5 reenters fiber optic 377 and is returned to fiber junction 375. The light having Interacted with the blood is 
returned to the spectrometer by fiber optic 379. The returning light is then separated by wavelength by 
grating 319 and reflected onto array detector 321, capable of detecting frequencies of light in the visible 
and infrared frequency region. 

The light intensities at various wavelengths are then analyzed by processing unit 305 employing a 

10 appropriate algorithm. Instrument 301 will display the current values of pH, PCO2, [HCO3"], PO2 and O2 
sat. as well as the past history of said analyte in said displays 351-359. If a given analyte determination 
indicates that the analysis might be unreliable the outlier lights 361-369 will indicate such to the operator. 

Whereas the drawings and accompanying description have shown and described the preferred 
embodiment of the present invention, it should be apparent to those skilled in the art that various changes 

75 may be made in the form of the invention without affecting the scope thereof. 

Claims 

1, A method of determining non-invasively and in vivo at least two of the five blood gas parameters (i.e., 
20 pH, [HCO3"], PCO2, PO2. and O2 sat.) in a human, said method comprising steps of: 

a. generating three or more different wavelengths of light, said wavelengths being in the range of 
500 nm to 2,400 nm; 

b. irradiating In vivo and non-invasively blood containing tissue with said wavelengths of light, so that 
there is differential attenuation of at least some intensities of said wavelengths by said blood 

25 containing tissue as a function of said wavelengths, said differential attenuation causing intensity 

variations of said wavelengths incident from said blood containing tissue as a function of said 
wavelengths, said tissue and said unknown levels of PH, [HCO3"], PCO2, PO2, and O2 sat.; 

c. measuring said Intensity variations from said blood containing tissue to obtain a set of intensity 
variations v. wavelengths; 

30 d. generating at least one calibration model, said model being a function of conventionally measured 

blood gas parameters obtained from a set of calibration samples and intensity variations v. 
wavelengths data obtained from irradiating said set of calibration samples with 3 or more different 
wavelengths of light, said wavelengths being in the range of 500 nm to 2,400 nm, said wavelengths 
being the same as used to irradiate said tissue; and 

35 e. calculating the values of said unknown blood gas parameters in said blood containing tissue from 

said measured intensity variations from said tissue utilizing at least one algorithm and said at least 
one calibration model, said calculation using 2 or more variables. 

2. The method as set forth in claim 1, further including the step of determining whether said intensity 
40 variations v. wavelengths from said tissue having said unknown blood gas parameters is an outlier, said 

determination including the step of comparing said Intensity variations v. wavelengths response from 
said tissue to said Intensity variations v. wavelengths from said model, to determine a measure of the 
magnitude of the difference between said tissue intensity variations v. wavelengths and said model 
intensity variations v. wavelengths. 

45 

a The method as set forth In claim 1, further including the step of pretreatment of said measured intensity 
variations and said measured blood gas parameters from said tissue, said pretreatment selected from 
the group including trimming, frequency, selection, centering, scaling, normalizing, taking first or higher 
order derivatives, smoothing, Fourier transfomning, linearization and transformation. 

50 

4. The method as set forth In claim 1. wherein said step of calculating is performed by a multivariate 
algorithm using 2 or more variables. 

5. The method as set forth in claim 4, wherein said algorithm is selected from the group including PLS, 
65 PLS2, PCR, CLS. Q-matrix, ridge regression, cross-correlation, Kalman filtering, MLR. neural networks, 

and continuum methods. 
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6. The method as set forth in claim 1 , wherein said calculating of said values of said unknown blood gas 
parameters is accomplished by utilizing at least two algorithms and wherein at least one of said two 
algorithms is a multivariate algorithm utilizing 2 or more variables. 

5 7. The method as set forth In claim 1, wherein said calibration model includes a calibration model for each 
blood gas parameter being determined by use of said differential attenuation of at least some 

intensities of said wavelengths. 

8. The method as set forth in claim 1, wherein said blood gas parameters are determined for arterial 
10 blood. 

9. The method as set forth in claim 8, wherein said step of irradiating said blood containing tissue 
includes the steps of irradiating during the diastolic portion of the cardiac cycle of said human and 
Inradiating said blood containing tissue during the systolic portion of said cardiac cycle, and wherein 

15 said step of measuring includes the step of measuring said intensity variations during said diastolic 
portion to obtain a diastolic set of intensity variations v. wavelengths, and the step of measuring said 
intensity variations during said systolic portion to obtain a systolic set of intensity variations v. 
wavelengths, and further including the step of pretreating said diastolic and systolic intensity variations 
V. wavelengths to determine the measure of change between said diastolic and systolic sets, and 

20 calculating said unknown blood gas parameters with the use of said measure of change. 

10. A method of determining invasively and in vivo at least two of the five blood gas parameters (i.e., pH, 
[HCO3"]. PCO2. PO2, and O2 sat.) in a human, said method comprising steps of: 

a. generating three or more different wavelengths of light, said wavelengths being in the visible and 
25 Infrared region; 

b. irradiating arterial blood with said wavelengths of light, so that there is differential attenuation of at 
least some intensities of said wavelengths by said blood as a function of said wavelengths, said 
differential attenuation causing intensity variations of said wavelengths incident from said blood as a 
function of said wavelengths, and said unknown levels of pH, [HCO3""], PCO2. PO2, and O2 sat. in 

30 said blood; 

c. measuring said intensity variations from said blood to obtain a set of Intensity variations v. 
wavelengths; 

d. generating at least one calibration model, said model being a function of conventionally measured 
blood gas parameters obtained from a set of calibration samples and intensity variations v. 

35 wavelengths data obtained from irradiating said set of calibration samples with 2 or more different 

wavelengths of light, said wavelengths being the same as used to irradiate said blood; and 

e. calculating the values of said unknown blood gas parameters in said blood from said measured 
intensity variations from said blood utilizing at least one algorithm and said at least one calibration 
model, said calculation using 2 or more variables. 

40 

11. A method of determining non-invasively and in vivo at least one of the five blood gas parameters (pH, 
[HCOa"], PCO2, and PO2) in a human, said method comprising steps of: 

a. generating two or more different wavelengths of light, said wavelengths being in the range of 500 
nm to 2,400 nm; 

45 b. irradiating in vivo and non-invasively blood containing tissue with said wavelengths of light, so that 

there is differential attenuation of at least some Intensities of said wavelengths by said blood 
containing tissue as a function of said wavelengths, said differential attenuation causing intensity 
variations of said wavelengths incident from said blood containing tissue as a function of said 
wavelengths, said tissue and said unknown levels of pH, [HCO3"], PCO2, and PO2; 

50 c. measuring said intensity variations from said blood containing tissue to obtain a set of intensity 

variations v. wavelengths; 

d. generating at least one calibration model, said model being a function of a conventionally 
measured blood gas parameter obtained from a set of calibration samples and intensity variations v. 
wavelengths data obtained from irradiating said set of calibration samples with two or more different 

55 wavelengths of light, said wavelengths being in the range of 500 nm to 2,400 nm. said wavelengths 

being the same as used to Inradiate said tissue; and 

e. calculating the value of said unknown blood gas parameter in said blood containing tissue from 
said measured intensity variations from said tissue utilizing at least one algorithm and said at least 
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one calibration model, said calculation using 1 or more variables. 

12. A non-Invasive and in vivo blood gas monitor for accurately and precisely determining at least two 
blood gas parameters in a human, said monitor comprising: 

5 a. means for generating a plurality of different wavelengths of light, said wavelengths being in the 

range of 500 nm to 2,400 nm; 

b. means for directing at least a portion of said wavelengths of light to a section of blood containing 
tissue of a human having unknown blood gas parameters, said means for directing including means 
for coupling to said human; 

10 c. means for collecting at least a portion of said wavelengths of light which are directed from said 

blood containing tissue to said means for collecting, said means for collecting Including some 
coupling means; 

d. means, connected to said means for collecting, for measuring the Intensity of at least two 
wavelengths collected by said means for collecting; 
75 e. means, coupled to said means for measuring, for storing said measured intensity variations v. 

wavelengths; 

f. means for storing at least one calibration model, said model being generated by an algorithm 
using 1 or more variables, said model being a function of said known blood gas parameters from 
said set of calibration samples and the intensity variations v. wavelengths obtained from irradiating 

20 and calibration samples with wavelengths of light in the range of 500 nm to 2,400 nm; 

g. means for storing said algorithm which utilizes said calibration model and said stored intensity 
variations v. wavelengths for determination of said unknown blood gas parameters; 

h. a microprocessor, coupled to said means for storing said calibration model, said means for storing 
said algorithm and said means for storing said intensity variations v. wavelengths, for calculating said 

25 unknown blood gas parameters; and 

i. means for Indicating said determined blood gas parameters. 

13. An invasive, in vivo blood gas monitor for accurately and precisely determining at least two blood gas 

parameters in a human said monitor comprising: 
30 a. means for generating a plurality of different wavelengths of light, said wavelengths being in the 

visible and infrared regions of the spectrum; 

b. means for directing at least a portion of said wavelengths of light into the blood of a human having 
unknown blood gas parameters, said means for directing including means for coupling to said 

human; 

35 c. means for collecting at least a portion of said wavelengths of light which are directed from said 

blood to said means for collecting, said means for collecting including said coupling means; 

d. means, connected to said means for collecting, for measuring the intensity of at least two 
wavelengths collected by said means for collecting; 

e. means, coupled to said means for measuring, for storing said measured intensity variations v. 
40 wavelengths; 

f. means for storing at least one calibration model, said model being generated by a algorithm using 
1 or more variables, said model being a function of said known blood gas parameters from said set 
of calibration samples and the intensity variations v. wavelengths obtained from irradiating said 
calibration samples with wavelengths of light in visible and infrared region, 

45 g. means for storing said algorithm which utilizes said calibration model and said stored intensity 

variations v. wavelengths for determination of said unknown blood gas parameters; 
h. a microprocessor, coupled to said means for storing said calibration model, said means for storing 
said algorithm and said means for storing said intensity variations v. wavelengths, for calculating said 
unknown blood gas parameters; and 

50 I. means for indicating said calculated blood gas parameter. 
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